Direct Observation of Dynamical Switching between Two Driven Oscillation States of 

a Josephson Junction. 
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We performed a novel phase sensitive microwave reflection experiment which directly probes the 
dynamics of the Josephson plasma resonance in both the linear and non-linear regime. When the 
junction was driven below the plasma frequency into the non-linear regime, we observed for the first 
time the transition between two different dynamical states predicted for non-linear systems. In our 
experiment, this transition appears as an abrupt change in the reflected signal phase at a critical 
excitation power. 



As first understood by Josephson, a superconduct- 
ing tunnel junction can be viewed as a non-linear elec- 
trodynamic oscillator 0. The tunneling of Cooper 
pairs manifests itself as a non-linear inductance that 
shunts the linear junction self-capacitance Cj, formed 
by the junction electrodes and the tunnel oxide layer. 
The constitutive relation of the non-linear inductor 
can be written as I(t) — IosmS(t), where I(t), 
5 (t) = J_ dt'V(t')/ipo and V(t) are the current, gauge- 
invariant phase-difference and voltage corresponding to 
the inductor, respectively, while the parameter Jo is the 
junction critical current. Here tpo — K/2e is the reduced 
flux quantum. For small oscillation amplitude, the fre- 
quency of oscillation is given for zero bias current by 
the so-called plasma frequency up = 1/^/LjCj where 
L.7 = Po/Io is the effective junction inductance. As the 
oscillation amplitude increases, the oscillation frequency 
decreases, an effect which has been measured in both the 
classical and quantum regime 0, IE 0, IE IE • However, 
a more dramatic non-linear effect should manifest itself 
if the junction is driven with an AC current isinwt at 
a frequency u> slightly below cup. If the quality factor 
Q = Cjujp/Re[Z~ 1 (u;p)] is greater than y/3/2a, where 
Z(ojp) is impedance of the junction electrodynamic envi- 
ronment and a = 1—uj/lup the detuning parameter, then 
the junction should switch from one dynamical oscillation 
state to another when i is ramped above a critical value 
id lEl- For i < i c i, the oscillation state would be low- 
amplitude and phase-lagging while for i > id, the oscil- 
lation state would be high-amplitude and phase-leading. 
This generic non-linear phenomenon, which we refer to as 
"dynamical switching" , is reminiscent of the usual "static 
switching" of the junction from the zero-voltage state 
to the voltage state when the DC current bias exceeds 
the critical current Iq 0]. However, an important dis- 
tinction between dynamical and static switching is that 
in dynamical switching, the phase particle remains con- 
fined to only one well of the junction cosine potential 
U(S) = —ipolocos(5), and the time-average value of S is 
always zero. The junction never switches to the "nor- 
mal" state, and thus no DC voltage is generated. Also, 



for dynamical switching, the current id depends both 
on Q and on the detuning a. In this Letter we report 
the first direct observation of this dynamical switching 
effect in a Josephson junction. The parallel with static 
switching suggests that dynamical switching can be used 
for amplification with the added advantage that no en- 
ergy dissipation occurs in the junction chip, a desirable 
feature for the readout of superconducting qubits 0] . 

Typical junction fabrication parameters limit the 
plasma frequency to the 20 - 100 GHz range where tech- 
niques for addressing junction dynamics are inconvenient. 
We have chosen to shunt the junction by a capacitive ad- 
mittance to lower the plasma frequency by more than an 
order of magnitude and attain a frequency in 1-2 GHz 
range (microwave L-band). In this frequency range, a 
simple on-chip electrodynamic environment with mini- 
mum parasitic elements can be implemented, and the 
hardware for precise signal generation and processing is 
readily available. In our experiment, we directly measure 
the plasma resonance in a coherent microwave reflection 
measurement. Unlike previous experiments which mea- 
sured only the microwave power absorption at the plasma 
resonance IE > we a l so measure the phase cj) of the 
reflected microwave signal. Thus, we can detect the char- 
acteristic signature of the transition between different os- 
cillating states of the junction - a change of oscillation 
phase relative to the drive. Note that the phase 4> of the 
reflected signal which probes the phase of the oscillation 
state should not be confused with the junction gauge- 
invariant phase difference 8. 

In the first step of sample fabrication, a metallic un- 
derlayer - either a normal metal (Au, Cu) or a super- 
conductor (Nb) - was deposited on a silicon substrate to 
form one plate of the shunting capacitor, followed by the 
deposition of an insulating SiaN4 layer. Using e-beam 
lithography and double-angle shadow mask evaporation, 
we subsequently fabricated the top capacitor plates along 
with a micron sized AI/AI2O3/AI tunnel junction. The 
critical current of the junction was in the range Iq = 
1 — 2 /iA. By varying both the dielectric layer thickness 
and the pad area, the capacitance C was varied between 



2 




FIG. 1: Schematic of the measurement setup. Thick lines 
correspond to 50 Q coaxial transmission lines. A lumped el- 
ement model for the junction chip and measurement line is 
shown. The two ideal current sources actually represent ex- 
ternal sources. 



16 and 40 pF. Sample parameters are listed in Table I. 

The junction + capacitor chip is placed on a microwave 
circuit-board and is wire-bonded to end of a coplanar 
stripline which is soldered to a coaxial launcher affixed 
to the side wall of the copper sample box. We anchor the 
RF leak-tight sample box to the cold stage of a 3 He re- 
frigerator with base temperature T — 280 mK. The mea- 
surement setup is schematically shown in Figure 1. Mi- 
crowave excitation signals are generated by a HP 8722D 
vector network analyzer and coupled to the sample via 
the -13 dB side port of a directional coupler. The re- 
flected microwave signal passes through the direct port 
of the coupler, and is amplified first using a cryogenic 
1.20 — 1.85 GHz HEMT amplifier with noise temperature 
TV = 4K before returning to the network analyzer. A 
DC bias current can be applied to the junction by way 
of a bias tee. We use cryogenic attenuators and isola- 
tors on the microwave lines in addition to copper-powder 
and other passive filters [I| on the DC lines to shield the 
junction from spurious electromagnetic noise. 

We locate the linear plasma resonance by sweeping 
the excitation frequency from 1 to 2 GHz and measur- 
ing the reflection coefficient T(oj) — i^/ii e?^ — (Z(uj) — 
Zq)/(Z(lu) + Zo), where Zq — 50 f2 is the character- 
istic impedance of our transmission lines and Z(lo) is 
the impedance presented to the analyzer by the chip 
and the measurement lines. For an ideal LC resonator 
without intrinsic dissipation, we expect a phase shift 
= w > Wp - 0w<cw P = 27r, which we verified by placing 
a chip capacitor and an inductive wire bond in place of 
the junction chip. An important aspect of our experiment 
is that Q is now determined by the ratio Zq/Zj ~ 10, 
where Zj — WLj/JCj + C) and not by the intrinsic 
junction losses which are negligible. An excitation power 
P = i 2 Z /4 « — 120dBm (1 fW) corresponding to a cur- 



Sample 


Lj(nH) 


oj p /27r(GHz) 


C(pF) 


Ls(nH) 


1 


0.28 


1.18 


39±1 


0.20±.02 


2 


0.18 


1.25 


30±4 


0.34±.04 


2a 


0.17 


1.66 


18±1 


0.32±.02 


3 


0.32 


1.64 


16±1 


0.27±.02 


4 


0.38 


1.81 


19±1 


0.026±.02 


5 


0.40 


1.54 


19±1 


0.15±.02 



TABLE I: Sample parameters. Lj — tpo/Io and u> p are mea- 
sured values. C and Ls are fit values to the data. Samples 
1,2 and 2a have a 100 nm thick Au underlayer, sample 3 has 
a 50 nm thick Nb underlayer, sample 4 has a 1 /im thick Cu 
underlayer, and sample 5 has a 200 nm thick Nb underlayer. 



rent i — 9 nA <C Iq keeps the junction in the linear 
regime. 

In Figure 2, we present the reflected signal phase 4> as a 
function of excitation frequency for sample 5. In order to 
remove the linear phase evolution associated with the fi- 
nite length of the measurement lines, we have subtracted 
from our measurement in the superconducting state the 
reflection coefficient measured with the junction in the 
normal state. The point where <j) = is the linear-regime 
plasma frequency. For sample 5, uj p / 2ir — 1.54 GHz. 

The precise frequency and critical current dependence 
of the reflected signal phase of our samples can be ac- 
counted for by a 3-element model for the electrodynamic 
environment seen by the junction. This lumped ele- 
ment model is shown in the lower right corner of Fig- 
ure 1. The parasitic inductance Ls and resistance Rs 
model the non-ideality of the shunting capacitor C. They 
arise from the imperfect screening of currents flowing in 
the capacitor plates and the finite conductivity of these 
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FIG. 2: Normalized reflected signal phase <f> as a function 
of excitation frequency for sample 5. The open circles are 
measured data for L j — 0.40 nH. The solid line is calculated 
from the equivalent circuit model shown in the inset. The 
magnitude of the reflected signal is unity. 
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plates. The plasma frequency in the linear regime is de- 
termined by the total inductance Lj + Lg and capaci- 
tance C e f f = Cj + C ~ C, and is given by the following 
relation: 

(-) 2 = C(L ; + L s ) = ^f + CL S . 

We thus plot 1/ojp versus 1/Iq = Lj/cpo in Figure 3 for 
samples 1, 2, 4 and 5. As the critical current is decreased 
by applying a magnetic field, the junction inductance in- 
creases, and the plasma frequency is reduced. For each 
sample, a linear fit to the data yields the values of C and 
Ls- The fit values for C agree well with simple estimates 
made from the sample geometry. Samples 1 and 2 have 
nominally the same capacitance but different critical cur- 
rent, and hence lie approximately on the same line in Fig- 
ure 3. A total of four capacitive pads were used to make 
the shunting capacitor in samples 1 and 2, and after ini- 
tial measurements, we scratched off two of the pads from 
sample 2 to obtain sample 2a. The fit parameters for 
sample 2a indicate that the capacitance is indeed halved. 
For samples with a thin underlayer (1,2 and 3), a stray 
inductance in the range Ls = 0.20 — 0.34 nH is observed. 
For samples 4 and 5 with a significantly thicker under- 
layer, Ls was reduced to 0.026 nH and 0.15 nH respec- 
tively. This behavior is consistent with the calculated 
screening properties of our thin films. To verify that the 
values of C and Ls were not affected by the magnetic 
field used to vary Iq, we varied Lj by applying a DC 
bias current Q at zero magnetic field. The resonance 
data obtained by this method agrees with the magnetic 
field data. With a single set of parameters, we can ac- 
curately fit both the position of the resonance and its 
lineshape. Additionally, we find Rs — 0.8 il for samples 
1 and 2, R s = 0.02 for sample 4, and R s — for the 
superconducting samples 3 and 5. The fit for sample 5 is 
shown as a solid line in Figure 2. Finally, we have inde- 
pendently verified the effect of the shunting capacitor on 
the plasma resonance by performing resonant activation 
experiments 

In order to study the non-linear regime of the plasma 
resonance, we measured the reflection coefficient as a 
function of frequency for increasing power. We present 
the data for sample 5 as a two dimensional color plot in 
the lower panel of Figure 4 in which each row is a sin- 
gle frequency sweep, similar to Figure 2. For small ex- 
citation power, we recover the linear plasma resonance 
at 1.54 GHz, shown as a yellow line corresponding to 
= 0. As the power is increased above — 115 dBm, the 
plasma frequency decreases, as is expected for a junc- 
tion driven with large amplitude Q. The boundary be- 
tween the leading-phase region (green) and the lagging- 
phase region (red) therefore curves for high powers. This 
curvature has an interesting consequence: When we in- 
crease the drive power at a constant frequency slightly 
below the plasma frequency, the phase as a function of 
power undergoes an abrupt step, as predicted 8]. For 
yet greater powers (> —90 dBm), we encounter a new 
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FIG. 3: Inverse square of the plasma frequency 1/Wp as a func- 
tion of the inverse critical current l/Io = Lj/tpo for samples 
1,2,4 and 5. Solid lines are linear fits to the data correspond- 
ing to the model of Fig. 1, with a single best fit line drawn 
for samples 1 and 2 which nominally differ only in Jo- 

dynamical regime (black region in Figure 4) where 5 
appears to diffuse between the wells of the cosine po- 
tential. When measuring the AC junction resistance at 
/ = 221 Hz in presence of the microwave drive, a finite 
resistance was only observed in the black region. To ver- 
ify that this phenomenon was not due to junction heat- 
ing, we have varied the ramp time of the power and no 
change in the transition power was observed. In the lower 
panel of Figure 4, we illustrate the sequence of dynam- 
ical transitions by plotting as a function of incident 
power at uj/2tt = 1.375 GHz. For P < -102 dBm, the 
phase is independent of power (5 oscillates in a single 
well in the harmonic- like, phase-leading state, letter A). 
For -102 dBm < P < -90 dBm, the phase evolves with 
power and S still remains within the same well, but os- 
cillates in the anharmonic phase-lagging state (letter B). 
Finally, for P > —90 dBm, the average phase of the re- 
flected signal saturates to -180 degrees, corresponding to 
a capacitive short circuit (letter C). This last value is 
expected if 6 hops randomly between wells, the effect 
of which is to neutralize the Josephson inductive admit- 
tance. 

The value of the switching current i c \ for the A-B tran- 
sition, which is a function of both the detuning a and 
power P, is in good qualitative agreement with the an- 
alytical theory which retains only the first anharmonic 
term in the cosine potential |8J. For instance, the slope 
of the A-B transition line at the linecut in Figure 4, 
dP (dBm) / da(%) — 0.8 for the experiment while we cal- 
culate its value to be 0.7. Furthermore, in measurements 
in which the power is ramped in less than 100 ns, we ver- 
ified that the transition between dynamical states is hys- 
teretic, another prediction of the theory. These results 
will be presented in a later publication. To explain the 
complete frequency and power dependence of the transi- 
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FIG. 4: Normalized reflected signal phase (j> (wrap-around 
color scheme) as a function of excitation frequency u/2n and 
excitation power P. In the lower panel, a vertical slice taken 
at ui/2n = 1.375 GHz (dashed line) shows the abrupt transi- 
tion between two oscillation states of the system. The upper 
panel is the result of numerical simulations. 



tions shown in the lower panel of Figure 4, we have per- 
formed numerical simulations by solving the full circuit 
model of the lower corner of Figure 1, including the ex- 
act junction non-linear constitutive relation. The result 
of this calculation is shown in the upper panel of Figure 
4. It correctly predicts the variation of the plasma fre- 
quency with excitation power, and the boundaries of the 
phase diffusion region. The agreement between theory 
and experiment is remarkable in view of the simplicity 
of the model used with no adjustable parameters, and 
only small differences in exact shape of region bound- 
aries are observed. It is important to mention that the 
overall topology of Figure 4 is unaffected by changes in 
the parameter values within the bounds of Table I. 

In conclusion, we have performed a novel, phase- 
sensitive, microwave experiment demonstrating that the 
Josephson plasma oscillation can switch between the two 
dynamical states predicted for a driven non-linear sys- 
tem. At the critical excitation power we observe an 
abrupt change in the reflected signal phase. In accor- 
dance with the analytical theory and numerical simu- 
lations, we have observed that this critical power is a 
strong function of the junction critical current. This 
phenomenon can therefore be applied to the detection 
of small relative variations in critical current in the same 
manner that the usual DC SQUID, biased in the vicinity 
of the critical current, transforms a flux-induced variation 
in critical current into a DC voltage 1 1 1| . Furthermore, 
following the methodology invented for experiments with 
trapped electron systems [12|, we can use this dynami- 
cal switching as the basis for a single-shot and latching 
qubit readout. Since the measurement of phase is purely 
dispersive, this new readout would have advantage the of 
eliminating dissipation at the chip level, a limiting factor 
for superconducting qubits. 
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